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ABSTRACT 


The  purpose  of  this  study  was  to  gain  a  better  understanding  of  the 
parameters  affecting  crack  growth  direction  in  unidirectional  composite 
materials.  To  achieve  this,  the  effect  of  anisotropy  and  biaxial  load¬ 
ing  on  the  direction  of  crack  growth  in  unidirectional  off-axis  compos¬ 
ite  materials  were  investigated.  Specific  emphasis  was  placed  on  defin¬ 
ing  the  crack  tip  stress  field  and  finding  a  consistent  criterion  for 
predicting  the  direction  of  crack  growth. 

Two  models  are  presented  to  predict  the  crack  tip  stress  field,  an 
anisotropic  elasticity  solution  and  a  singular  isoparametric  finite 
element  formulation.  After  defining  the  crack  tip  stress  field,  three 
crack  extension  direction  criteria,  the  normal  stress  ratio,  the  tensor 
polynomial  and  the  strain  energy  density  criterion,  were  applied  to 
predict  the  direction  of  crack  extension. 

The  theoretically  predicted  crack  extension  directions  were  then 
compared  with  experimental  results.  After  comparison,  it  was  determined 
that  only  the  normal  stress  ratio  criterion  correctly  predicts  the 
direction  of  crack  extention. 
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Chapter  1 
INTRODUCTION 

In  the  last  twenty  years,  a  desire  for  structural  materials  with 
greater  strength  to  weight  ratios  has  led  to  the  development  of  fiber 
reinforced  composites.  Early  on,  composites  were  used  in  applications 
that  were  not  structurally  critical.  For  those  applications,  character¬ 
ization  of  the  mechanical  properties  of  the  material  was  often  all  that 
was  necessary  before  application.  Now,  when  composites  are  used  for 
main  load  carrying  components  in  structures,  an  understanding  of  the 
failure  modes  of  the  composite  is  necessary.  Of  particular  interest  is 
the  fracture  response  of  composites,  i.e.,  the  influence  of  initial 
flaws  and  flaws  which  develop  during  service,  on  the  structural  integ¬ 
rity  of  the  component. 

1.1  Motivation  and  Statement  of  Purpose 

A  fundamental  problem  in  predicting  the  failure  of  laminated  com¬ 
posite  materials  is  an  understanding  of  the  direction  of  crack  growth  in 
the  individual  laminae.  The  importance  of  the  direction  of  crack  growth 
on  the  failure  mode  of  the  laminate,  is  shown  in  Fig.  1.0  [!]• 

Clustered  [{0)2/(-e)2]s  graphite-epoxy  laminates  failed  in  a  pure  matrix 
mode  (delamination  and  either  intralaminar  matrix  cracking  or  fiber 
matrix  debonding).  In  contrast,  the  alternating  [(■*'0/"0)2]5  laminates 
exhibited  fiber  breakage  in  half  of  the  plies  and  either  matrix  cracking 
or  fiber  matrix  debonding  in  the  others;  there  was  no  delamination  in 
the  alternating  laminates.  The  mode  of  failure  has  a  significant  effect 
on  the  strength  of  the  laminate.  The  strength  of  the  10®  and  30®  alter- 
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nating  laminates  was  30  and  50  percent  greater,  respectively,  than  the 
strength  of  the  clustered  laminates  Cll»  Hence,  understanding  the 
parameters  that  affect  the  laminate  failure  mode,  particularly  those 
Influencing  the  direction  of  crack  growth  in  the  lamina,  is  of  critical 
importance  in  predicting  the  fracture  response  of  laminates. 

The  purpose  of  this  study  is  to  gain  a  better  understanding  of  the 
parameters  affecting  crack  growth  direction  in  the  lamina.  To  achieve 
this,  the  effect  of  anisotropy  and  biaxial  loading  on  the  direction  of 
crack  growth  in  unidirectional  off-axis  composite  materials  will  be 
investigated.  Specific  emphasis  will  be  placed  on  defining  the  crack 
tip  stress  field  and  finding  a  consistent  criterion  for  predicting  the 
direction  of  crack  growth. 

1.2  Literature  Review 

Several  researchers  have  addressed  fracture  of  unidirectional 
composites  in  the  past.  Wu  [2]  and  Sih  et  a1 .  [3]  have  Independently 
investigated  the  fracture  response  of  unidirectional  glass/epoxy  com¬ 
posites,  They  have  shown  success  in  predicting  the  critical  load  for 
mode  1,  mode  2,  and  mixed  mode  problems,  with  the  restriction  that  the 
initial  crack  be  aligned  in  the  fiber  direction.  Additional  work  by  Wu, 
some  of  which  applies  to  the  fracture  of  composite  laminates  can  be 
found  in  [4-6],  while  Sih  summarizes  much  of  the  current  research 
efforts  in  composite  fracture  in  [7],  The  Wu  and  Sih  theories  will  be 
discussed  more  fully  later  in  this  text. 

More  recently,  Goree  and  Jones  [8]  examined  the  fracture  behavior 
of  unidirectional  boron/aluminum  composites.  In  their  analysis,  Goree 


and  Jones  treat  the  material  as  an  Inhomoyeneous  anisotropic  body  con¬ 
taining  broken  fibers,  undamaged  fibers  and  matrix  material.  They 
assume  that  the  fibers  carry  all  of  the  axial  load  and  that  the  matrix 
transfers  load  from  a  fiber  break  to  an  adjacent  unbroken  fiber  by 
shear.  The  Initial  crack  is  modeled  by  an  arbitrary  number  of  fiber 
breaks  and  the  direction  of  crack  extension  is  predicted  by  a  maximum 
stress  criterion. 

These  researchers  have  made  significant  contributions  to  under¬ 
standing  the  fracture  response  of  unidirectional  composites.  However,  a 
more  general  analysis  is  needed.  In  the  preceding  models,  numerous 
assumptions  have  been  made  on  the  orientation  of  the  initial  crack,  the 
scope  of  the  analysis  (microscopic  versus  macroscopic),  and  on  the  final 
mode  of  failure.  In  order  to  provide  a  more  general  theory,  Buczek  and 
Herakovich  [9,10]  studied  the  direction  of  crack  growth  and  the  associ¬ 
ated  energy  release  rates  in  a  unidirectional  off-axis  lamina  and  edge 
delamination  in  a  [O/gojj  laminate.  They  treated  the  lamina  as  a  homo¬ 
geneous  anisotropic  material.  Using  the  finite  element  method  to  model 
the  crack  tip  stress  field  in  conjunction  with  various  crack  growth 
criteria,  Buczek  and  Herakovich  predicted  the  direction  of  crack  exten¬ 
sion  in  unidirectional  composite  tensile  coupons.  Though  the  analysis 
was  limited  to  tensile  coupons  subjected  to  uniaxial  displacement  load¬ 
ing,  no  restrictions  were  made  on  the  orientation  of  the  Initial 
crack.  Buczek  and  Herakovich  also  noted  several  inconsistencies  among 
the  crack  growth  direction  criteria  previously  presented.  As  a  result, 
they  have  hypothesized  a  model  to  predict  the  direction  of  crack  exten¬ 
sion  in  unidirectional  composites.  For  the  limited  number  of  cases 
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analyzed,  the  direction  of  crack  yrowth  predicted  by  their  new  model 
correlates  well  with  experimental  evidence. 


Chapter  2 

THEORETICAL  CONSIDERATIONS 

2.1  Crack  Growth  Direction  Criteria 

In  this  section  three  crack  growth  direction  criteria  are  presen¬ 
ted.  Though  these  criteria  can  also  be  used  to  predict  load  at  failure, 
the  main  emphasis  in  this  study  is  on  predicting  crack  extension  direc¬ 
tion.  The  criteria  are  applied  treating  the  unidirectional  composite 
lamina  as  a  homogeneous  anisotropic  material.  Salient  features  of  these 
criteria  must  be  the  applicability  to  mixed-mode  fracture  problems  and 
the  ability  to  account  for  the  anisotropic  nature  of  crack  growth  resis¬ 
tance  in  the  unidirectional  lamina  (i.e.,  accounting  for  the  fact  that 
it  is  more  difficult  to  break  fibers  than  matrix).  Three  crack  exten¬ 
sion  criteria  addressing  these  factors  have  been  presented  In  the  liter¬ 
ature.  They  are  the  Tensor  Polynomial  Failure  Criterion  [11],  the 
Strain  Energy  Density  Criterion  [12],  and  the  Normal  Stress  Ratio  Cri¬ 
terion  [10].  Though  the  Strain  Energy  Density  Criterion  does  not 
specifically  account  for  the  anisotropy  of  strength  in  composite  materi¬ 
als,  Sih  argues  in  [3]  that  the  criterion  is  applicable  at  least  in 
principle  to  anisotropic  fracture  problems. 

2.1.1  Tensor  Polynomial  Failure  Criterion 

The  Tensor  Polynomial  Criterion  Isa  phenomenological  failure 
theory  presented  by  Tsai  and  Wu  [11].  This  theory  is  based  on  the  exis¬ 
tence  of  a  failure  surface  in  stress  space  of  the  form: 

f(°i)  =  ''i  °i  ''ij  ®i  (2.1) 
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where  and  are  strength  tensors  of  second  and  fourth  order, 
and  0^  the  contracted  form  of  the  stress  tensor.  Expressions  for  F^  and 
F^j,  In  both  the  x-y  and  1-2  coordinate  systems  are  presented  in  Appen¬ 
dix  A. 

When  applying  the  Tensor  Polynomial  Criterion  to  fracture  problems, 
the  assumed  direction  of  crack  extension  is  the  radial  direction  of 
maximum  f(o^).  The  stress  components  are  those  determined  by  a 
continuum  mechanics  based  stress  analysis  and  must  be  evaluated  at  a 
finite  distance,  r^,  from  the  crack  tip.  The  crack  tip  coordinate 
system  is  shown  in  Fig.  2.0, 

A  graphical  application  of  the  Tensor  Polynomial  Criterion  to 
predict  the  direction  of  crack  growth  is  shown  by  Wu  in  [2].  This 
graphical  approach  is  presented  in  Fig.  2.1.  The  contour  of  the  stress 
vector  acting  on  radial  planes  emanating  from  the  crack  tip  is 
represented  by  S,  while  the  contour  of  the  strength  vector  on  the  same 
planes  is  represented  by  F.  The  predicted  direction  of  crack  extension 
is  labeled  <]>,.  It  should  be  noted  that  the  direction  of  does  not 
necessarily  correspond  to  the  direction  of  maximum  S  as  would  be  the 
case  for  isotropic  materials.  In  anisotropic  fracture  the  crack  exten¬ 
sion  resistance,  represented  by  F,  plays  an  important  role  in  determin¬ 
ing  4)^.  and  must  be  accounted  for. 

2.1.2  Strain  Energy  Density  Criterion 

The  Strain  Energy  Density  Criterion  is  phenomenological  in  nature 
and  is  based  on  variations  in  the  energy  stored  along  the  periphery  of  a 


where. 

F=  Contour  of  Strength  Vector 
S  =  Contour  of  Stress  Vector 

Pi’fidicted  Direction  of  Crack  Extension 

Fig.  2.1  Graphical  Application  of  the  Tensor  Polynomial  Criterion. 
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core  region  surrounding  the  crack.  Sih  presents  the  criterion  in  [12] 
for  isotropic  fracture  and  a  modified  form  for  application  to  anisotro¬ 
pic  fracture  in  [3]. 

Sih  defines  the  strain  energy  density  factor,  S,  as: 


aw  S 

3V  r 


(2.2) 


where  -jy  is  the  strain  energy  density  function  and  r  the  distance  from 
the  crack  tip.  Since  the  strain  energy  density  function  can  be 
expressed  in  terms  of  the  crack  tip  stresses  and  strains  as; 


aw  .  1 
W  7 


IVx 


o  c 

y  y 


+  T  Y  ) 

xy’xy' 


(2.3) 


an  expression  for  the  strain  energy  density  factor,, S,  can  be  obtained 
by  substitution.  This  expression  is: 


(2.4) 


The  fundamental  hypothesis  of  Sih  [3]  on  unstable  crack  growth  is: 

Crack  initiation  takes  place  in  the  radial  direction  corre¬ 
sponding  to  a  minimum  value  of  the  strain  energy  density  fac¬ 
tor,  i.e., 

2 

ly  »  0  and  -2-1  >  0  at  (j>  «  (2.5) 

3<{!> 

j 


n 


S1h  cautions  that  for  very  small  values  of  r  a  continuum  mechanics 
based  crack  tip  stress  analysis  is  invalid.  Hence  the  strain  energy 
factor  should  be  evaluated  at  a  finite  distance,  r^,  from  the  crack  tip, 
where  rg  is  of  the  same  order  of  magnitude  as  the  crack  tip  curvature. 
The  crack  tip  coordinate  system  for  this  criterion  is  similar  to  that  of 
the  Tensor  Polynomial  as  shown  in  Fig,  2.0. 

2.1.3  Normal  Stress  Ratio  Criterion 

Buczek  and  Herakovich  [lOJ  have  hypothesized  the  Normal  Stress 
Ratio  Criterion  as  a  phenomenological  crack  growth  direction  criter¬ 
ion.  The  model  assumes  that  the  direction  of  crack  extension  corre¬ 
sponds  to  the  direction  of  the  maximum  value  of  R(rQ,i)i),  defined  as 

(2.6) 

In  the  expression  for  R(rQ,iji),  corresponds  to  the  normal  stress 
acting  on  a  radial  plane  defined  by  at  a  given  distance  r^  from  the 
crack  tip,  and  T^^  is  the  tensile  strength  normal  to  the  ij>  plane. 

Since  the  tensile  strength  of  an  arbitrary  plane  is  difficult  if 
not  impossible  to  measure,  T^^  is  defined  in  a  consistent  manner  with 
the  tests  that  can  oe  performed.  To  meet  these  requirements,  a  mathe¬ 
matical  definition  of  T^^  must  satisfy  the  following  conditions: 

(1)  for  an  Isotropic  material,  T^^  must  be  independent  of 

(2)  for  crack  growth  parallel  to  the  fibers,  T,^  must  equal  the 
transverse  tensile  strength  Yj. 
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(3)  for  crack  growth  perpendicular  to  the  fibers,  must  equal 

99 

the  longitudinal  strength  Xi-. 

A  definition  satisfying  these  conditions  Is: 

V  “  ^  (2.7) 

where  0  Is  the  angle  from  the  plane  of  Interest  to  the  fiber  direc¬ 
tion.  These  parameters  along  with  the  crack  tip  coordinate  system  for 
the  Normal  Stress  Ratio  are  shown  In  Fig.  2.2. 

2.2  Crack  Tip  Stress  Field  Models 

Two  theoretical  models  to  predict  the  crack  tip  stress  field  In 
homogeneous  anisotropic  materials  are  presented  In  this  section.  The 
first  Is  an  anisotropic  elasticity  solution  for  a  blaxially  loaded 
center  cracked  Infinite  plate.  The  second  Is  a  singular  finite  element 
formulation  which  allows  for  variable  loading  and  geometries. 

2.2.1  Anisotropic  Elasticity  Solution 

The  stress  analysis  of  an  Infinite  homogeneous  anisotropic  plate 
containing  a  center  crack  can  be  directly  related  to  a  homogeneous 
anisotropic  plate  with  an  elliptic  hole.  By  reducing  the  minor  axis 
dimension  to  zero  and  evaluating  the  stress  potential  functions  In  the 
neighborhood  of  the  crack  tip,  Lekhnltskll’s  complex  variable  solution 
[13]  for  an  elliptic  hole  In  an  anisotropic  plate  can  be  applied  to 
anisotropic  fracture  problems.  Wu  presents  a  detailed  description  of 


this  application  In  [14]  along  with  equations  describing  the  crack  tip 
stresses  for  an  infinite  homogeneous  anisotropic  center  cracked  plate. 
Figure  2.3  shows  the  problem  under  consideration  along  with  the  crack 
tip  coordinate  system. 

The  governing  partial  differential  equation  for  the  problem  in 
terms  of  the  Airy's  stress  function  U  is: 

2A2g  3%  (2A,2  +  Agg)  2A.g  3^U  A,,  3% 

T  “  T  +  jj  y  - - y  + - ^  =0  (2,i 

3x  A22  ax  ay  >^zz 

where  A,j  are  components  of  the  compliance  tensor  for  plane  stress  or 
plane  strain  depending  on  the  analysis  desired. 

Assuming  U  =  e^  ^  the  characteristic  equation  for  (2.8)  takes 
the  form: 


AiiS  -  2Aj^gS  +  (2Aj2  Agg)S^  -  2A2gS  +  A22  =  0  (2.9) 

The  roots  of  the  characteristic  equation,  Sj  and  $2,  are  complex  and  are 
functions  of  the  material  properties  and  the  orientation  of  the  crack 
relative  to  the  material  principal  direction. 

The  solution  of  Eq,  2.8  can  be  obtained  in  terms  of  two  holomorphic 
functions,  and  *2(^2^’  following  form: 

U  =  z  Re{Fj(Zj)  +  F2(z2)|  (2.10) 


where 


dF,{z,) 


dF„(z5,) 
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and  the  complex  variables  are  defined  as: 

^1  “  X  +  Siy  .  Zg  “  X  +  S.^y  (2.12) 

The  stress  and  displacement  components  for  an  anisotropic  body 
expressed  in  terms  of  the  complex  potentials  ♦j(Zj)  and  #2(^2^ 

.  2Re  [s2*|(zj)  +  s2*’{z2)] 

Oy  »  2Re[*j(z^)  +  *2(22)]  (2.13) 

T^y  =  -2Re[Sj*|(Zj)  +  82*2(22)] 


u  »  2Re[pj*j(Zj)  +  P2*2(z2)] 

V  =  2Re[qj*j^{Zj)  +  q2*2(22)] 

where 

Pi  =  +  Aj2  -  ^16^1 

P2  “  Aj^s|  +  A^2  -  ^16^2 
„  _  ^2^1  ^22  ■  ^^26^1 

_  Ai2$2  •»•  A22  -  A2g$2 


(2.14) 


(2.15) 
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Savin  [15]  presents  complex  stress  potentials  satisfying  the  bound¬ 
ary  conditions  for  a  crack  in  an  infinite  plate  subjected  to  ten¬ 
sion  o"  and  shear  t*  at  infinity  as  indicated  in  Fig.  2,3.  These 
complex  potential  functions  are; 


where 


I  a ("  +  $20  ) 

“  2(Sj  -  $2) 

.  -a(T-  +  S/) 

^2^^2^  “  “  2(S^  -  $2)  '^2 


Cj  a  rcos^  +  S^rsin^i 


C2  *  rcoscji  +  S2rsini(i 


(2.16) 


(2.17) 


These  functions  were  obtained  through  a  conformal  mapping  of  the  crack 
into  a  unit  circle  and  satisfy  the  boundary  conditions: 


0  ”  T  “0 

y  xy 


-a  <  X  <  a 


(2.18) 


To  examine  the  crack  tip  stress  distribution,  these  stress  poten¬ 
tials  are  evaluated  as  the  c.  approach  1, 

J 

lim  ,  (t“  ✓!  +  Sjo*  ✓!) 

c  +1  *1(^1)  - - - 

'’I  2(Sj  -  S2)/2r(cos<|.  +  S2Sin«) 


lim  J,,  ,  .  h“ /a  +  Sjo" /S’) 

C  +1  ^o(^o)  ”  . 

^2  ^  ^  ^  2(Sj  -  S2)/5F(cosF^r3^rrn?X 


(2.19) 
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The  potentials  are  evaluated  as  the  approach  1,  because  the  conformal 

w 

mapping  transforms  the  region  surrounding  the  crack  to  the  Interior  of  a 
unit  circle  whose  radius  Is  represented  by  c..  The  transformation  of  a 

J 

crack  Into  a  unit  circle  Is  shown  In  Fly.  2.4. 

Substitution  of  (2.19)  Into  (2.13)  and  (2.14)  gives  the  stress  and 
displacement  distribution  near  the  crack  tip. 


0 


X 


a"  /a  S, Sj  S,  S,  t*  /a 

/2r  (Sj  -  Sg)  /W 


1 

(Sj  -  Sg) 


^  ^  1  T  /a  1  ,  , 

—  R.  I - [-r-TTi)* - I - 

(Sj  -  $2)  /2f  (Sj  -  Sg)  <<(2 


o  /a  S,Sp  1  1  t“  /a  1  S,  So 

'x, - I  -  (-17  -  -ri  I  * - I - 't  -  t)  ) 

/?r  (S,  -  Sj)  tY  *¥  ^  (h  -  Sj)  *7  *7 


u  .  a-  /5i7  Ke  Its  -  1  s^)[SiP;4^-  ) 


(2.20) 


V  .  /SF  R.  I(s^  ‘  S2Pl»(? ) 


^2.  n  3'2 


/2ar  Re  s^t W -  Vl  ]  1 


where 


lij  =  cos(|>  +  Sj^slnij) 


ij»2  =  coS(ti  +  S2Sln<|> 


Z -Plane 


f-Plane 


Conformal  Transformations 

i*ii2 


Fig.  2.4  Conforml  Happing  of  Crack  Into  a  Unit  Circle 
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The  crack  tip  stresses  exhibit  a  singularity  of  l//r  as  in  the 
isotropic  case.  However,  the  magnitude  of  the  stresses  are  not  simply  a 
function  of  the  stress  intensity  factors  as  In  the  isotropic  case;  the 
quantities  and  $2  affect  the  magnitude  of  the  stresses.  This  is 
an  important  difference  between  anisotropic  and  isotropic  fracture 
mechanics.  In  anisotropic  fracture  mechanics,  the  magnitude  of  the 
crack  tip  stresses  are  a  function  not  only  of  the  applied  load,  specimen 
geometry  and  crack  length,  but  also  the  material  properties  and  the 
orientation  of  the  crack  relative  to  the  principal  material  direction. 

2,2.2  Singular  Finite  Element  Formulation 

From  the  anisotropic  elasticity  solution,  it  is  evident  that  near 
the  crack  tip  there  exists  a  singularity  in  the  stress  field  of 
order  1//F,  A  finite  element  with  a  similar  singularity  Incorporated 
into  the  equations  describing  the  stress  field  over  the  domain  of  the 
element  would  be  advantageous.  Barsouro  presents  such  a  singular  element 
in  [16]  {the  idea  of  applying  singular  isoparametric  elements  to  frac¬ 
ture  problems  was  presented  independently  by  Barsoum  [16]  and  by 
Henshel I  and  Shaw  [18]).  The  advantage  of  Barsoum's  formulation  is  that 
a  special  crack  tip  element  is  not  necessary.  The  1//F  singularity  can 
be  generated  in  quadratic  isoparametric  elements  by  positioning  the  mid¬ 
side  nodes  near  the  crack  at  the  quarter  points.  The  technique  is 
depicted  in  Fig.  2.5.  (A  derivation  of  the  singularity  is  presented  in 
Appendix  B).  Freese  and  Tracey  point  out  in  [17]  that  this  singularity 
is  present  in  both  the  natural  isoparametric  triangle  and  the  quadratic 
triangular  element  generated  by  collapsing  one  side  of  an  eight  noded 


Conventional  Triangular 
Element 


Singular  Triangular 
Element 


Fig.  2.5  Conventional  and  Singular  Isoparametric  Triangular  Elements. 
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Isoparametric  quadrilateral  element.  However,  as  Barsoum  shows  In  [16], 
to  obtain  the  l//r  singularity  In  the  collapsed  quadrilateral  the 
overlaping  nodes  must  be  constrained  to  have  the  same  displacements.  An 
obvious  advantage  of  the  natural  Isoparametric  triangle  Is  that  multi¬ 
point  constraints  of  overlapping  nodes  Is  not  required.  Another  advan¬ 
tage  shown  by  Freese  et  al.  [17],  Is  that  the  natural  isoparametric 
triangle  allows  for  distorted  triangular  shapes  {e.g.,  pie-shaped  ele¬ 
ments)  while  the  collapsed  quadrilateral  requires  the  element  edges  to 
remain  straight. 

Though  Barsoum  [16]  only  presents  applications  of  the  singular 
Isoparametric  element  In  Isotropic  fracture  mechanics,  the  element  does 
have  applications  In  anisotropic  fracture.  Since  Freese  et  al .  [17] 
have  shown  that  the  natural  Isoparametric  element  can  represent  a  radial 
displacement  variation  of  the  form: 

-  A^  +  B^/r  +  C^r  (2.21) 

and  since  the  displacements  from  the  anisotropic  elasticity  solution 
take  the  form: 


“l  =  A^/r  (2.22) 

the  natural  isoparametric  element  has  as  many  applications  in  anisotro¬ 
pic  fracture  as  in  Isotropic  fracture. 

A  survey  of  alternative  /inite  element  methods  in  fracture  mechan¬ 
ics  Is  given  by  Hilton  and  Sih  in  [19].  The  emphasis  of  their  paper  is 
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on  calculation  of  stress  Intensity  factors  using  finite  elements;  how¬ 
ever,  several  state-of-the-art  finite  element  formulations  .•or  solution 
of  fracture  problems  are  presented. 


Chapter  3 

CRACK  TIP  STRESS  ANALYSIS 

As  previously  mentioned,  one  of  the  objectives  of  this  work  Is  the 
prediction  of  crack  tip  stresses  in  anisotropic  materials  without 
restrictions  as  to  the  orientation  the  crack  relative  to  the  principal 
material  direction  or  the  type  of  loading  applied.  In  this  chapter,  the 
application  of  the  anisotropic  elasticity  solution  and  the  singular 
Isoparametric  finite  element  formulation  to  the  general  fracture  problem 
presented  In  Fig.  3.U  Is  discussed.  Note  there  are  no  restrictions  as 
to  the  orientation  or  the  crack,  defined  by  the  angle  a,  or  that  of  the 
principal  material  (fiber)  direction  defined  by  3.  Also,  complete 
biaxial  loading  Is  allowed;  the  far  field  loads  can  all  be  specified 
Independently.  The  only  constraint  Is  that  the  crack  Is  assumed  to  have 
a  finite  height.  This  disallows  the  specification  of  far  field  loadings 
that  cause  crack  closure. 

3.1  Application  of  the  Anisotropic  Elasticity  Solution 

Valid  leadings  for  the  elasticity  solution  presented  In  Section 
(2.2.1)  are  shown  In  Fig.  2.3  and  consist  of  a  uniform  tensile 

e 

00  00 

stress  a  and  a  shear  stress  x  .  Though  the  solution  allows  for  Inde¬ 
pendent  specification  of  the  material  principal  direction  defined  by  the 
angle  0,  the  crack  Is  assumed  to  be  parallel  to  the  x-ax1s.  This  Is  not 
as  general  as  the  problem  defined  In  Fig.  3.0  which  allows  for  the 
specification  of  five  Independent  values  o  ,  o  ,  x  ,  o  and  6.  However, 
by  transforming  the  far  field  stresses  to  a  crack  tip  coordinate  system 
and  realization  that  the  far  field  stress  parallel  to  the  crack  does  not 
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contribute  to  the  singularity  (this  fact  is  discussed  in  Wu  [14])  the 
elasticity  solution  can  be  applied  to  problems  similar  to  that  of  Fiy. 

3.0. 

The  details  of  this  transformation  and  superposition  are  shown  in 
Fig.  3.1.  First,  the  far  field  stresses  of  Fig,  3.1a  are  transformed  to 
a  crack  tip  coordinate  system  defined  by  the  angle  a  shown  in  Fig.  3,1b 
(the  crack  tip  coordinate  system  is  represented  by  the  axes  5  and  n). 

In  the  t-n  coordinate  system,  the  effect  of  o*  and  on  the  crack  tip 
stress  field  can  be  evaluated  by  applying  Eq.  2,20.  The  stresses  are 
separated  into  singular  and  nun>singular  components  as  shown  in  Fig. 
3.1c.  Before  Eq,  2.20  can  be  applied,  the  components  of  the  compliance 
tensor,  A^j,  must  also  be  transformed  to  the  crack  tip  coordinate 
system.  Since  the  crack  was  originally  oriented  at  an  angle  a  and  the 
material  principal  direction  was  defined  by  the  angle  6,  the  compliance 
tensor  must  be  transformed  by  the  angle  (0-a),  Let  B  represent  the 
angle  (e-o). 

The  use  of  stress  transformation  along  with  the  fact  that  the  far 
field  stress  parallel  to  the  crack  dees  not  contribute  to  the  singular¬ 
ity,  allows  for  the  application  of  specific  elasticity  solutions  to 
general  fracture  problems.  In  the  preceding  problem,  an  elasticity 
solution  allowing  the  specification  of  two  independent  values  is  used  to 
solve  a  problem  requiring  the  specification  of  five  values. 


( 

1  ^ 


Fig.  3.1  Transformation  and  Superposition  of  Far  Field  Stresses. 
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3,2  Application  of  Sinyular  Isoparametric  Finite  Elements 

Since  many  fracture  problems  of  Interest  do  not  satisfy  the  Infi¬ 
nite  boundary  conditions  of  the  anisotropic  elasticity  solution,  an 
anisotropic  finite  element  fracture  mechanics  code  was  developed.  This 
code,  ANFRAC  (anisotropic  fracture  mechanics  code).  Incorporates  an 
Isoparametric,  quadratic,  triangular  element  and  the  Inherent  singular¬ 
ity  associated  with  this  element  when  the  mid-side  nodes  are  positioned 
at  quarter  points.  The  advantage  of  a  finite  element  model  is  the 
ability  to  analyze  specific  geometries  and  loading  conditions  that 
cannot  be  addressed  with  an  elasticity  solution, 

ANFRAC  is  a  linear,  two-dimensional,  quadratic.  Isoparametric 
finite  element  code.  The  outline  of  ANFRAC  is  similar  to  that  of  FEM2D, 
an  isoparametric  finite  element  code  presented  by  Reddy  [20],  A  des¬ 
cription  of  the  quadratic  Isoparametric  triangular  element  along  with 
practical  considerations  on  the  use  of  isoparametric  elements  is  given 
in  Bathe  [21], 

The  master  configuration  of  this  element  along  with  the  location  of 
the  Integration  stations  Is  shown  in  Fig,  3,2,  Numerical  integration 
over  triangular  domains  Is  discussed  by  Cook  [22],  The  Integration 
scheme  used  to  calculate  the  element  stiffness  matrix  Is  a  seven  point 
Gaussian  quadrature.  The  formula  for  the  numerical  Integration  is; 

7 

//  f  dA  =  I  w.  f{r,s)  (3,1) 

1=1 

where 

f  is  the  function  to  be  integrated  over  the  triangular  area  A, 
w^  is  the  weight  associated  with  the  1th  sampling  point. 
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and 

f{r,s)  is  the  value  of  f  evaluated  at  the  1th  sampling  point. 

The  interpolation  functions  and  weights  for  the  isoparametric  quadratic 
trangular  element  were  taken  from  Bathe  [21]  and  are  given  in  Table  3.0. 

As  discussed  in  Barsoum  [23],  the  ratio  of  element  length  to  half 
crack  length  (designated  i/a)  significantly  influences  the  crack  tip 
stress  field.  A  similar  effect  is  found  in  the  natural  Isoparametric 
element.  Freese  and  Tracey  [17]  present  stress  Intensity  factor  results 
for  an  isotropic  plate  with  an  edge  crack  using  the  natural  isoparamet¬ 
ric  trangular  element.  They  have  chosen  an  i/a  ratio  of  0.1,  and  pres¬ 
ent  data  that  is  independent  of  the  angle  41.  With  their  results  in 
mind,  a  quarter-symmetry  finite  element  mesh  was  developed  to  solve 
fracture  problems  involving  a  center  cracked  plate.  This  mesh  is  shown 
in  Fig.  3.3  and  an  exploded  view  of  the  crack  tip  elements  is  given  in 
Fig.  3.4.  The  t/a  ratio  of  this  mesh  is  0.1  and  the  ratic  of  2a/w  is 
also  0.1.  The  ratio  2a/w  is  a  measure  of  the  infinite  nature  of  the 
cracked  plate.  As  the  ratio  of  the  crack  length,  2a,  to  the  plate 
width,  w,  decreases,  the  infinite  boudary  conditions  of  the  elasticity 
solution  are  closer  to  being  satisfied.  A  2a/w  ratio  of  0.1  implies  the 
plate  width  is  an  order  of  magnitude  larger  than  the  crack  length.  For 
data  given  in  Hertzberg  [24],  an  isotropic  plate  of  this  size  adequately 
satisfies  the  infinite  boundary  conditions.  This  appears  to  be  true  for 
anisotropic  plates  based  on  data  provided  by  Smith  and  Kullinix  [25]. 

Since  the  mesh  shown  in  Fig.  3.3  is  a  quarter  symmetry  mesh  it  can 
only  be  used  for  problems  with  a  crack  orientation  of  0®  and  dis- 


Fig.  3.3 
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Fig. 3.4  for  Enlorgoment  of  Region  A-A 


Syitwetry  Mesh  of  Center  Cracked  Plate. 


Note:  See  Fig.  3.3  for  Relative  Location  of 
Region  A-A. 


Fig.  3.4  Exploded  View  of  Crack  Tip  Elements. 
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Table  3.0 

Interpolation  Functions  and  Numerical  Integration  Data 

Interpolation  Functions 

hj  “  2r^  -  3r  *  4rs  -  3s  +  2s^  +  1 
i\2  “  2r^  -  r 
hj  =  2s^  -  s 

»  -4r^  +  4r  -  4rs 
hg  =■  4rs 

hg  =  “4s^  +  4s  -  4rs 


Gauss  Point  Coordinates  and  Weights 


r  -  coordinates 

rj  »  0.10128  65073  235 

r2  =*  0.79742  69853  531 

''3  =  '“I 

r4  »  0.47014  20641  051 
••5  *  '“4 

rg  «  0.05971  58717  898 
r-f  =  0.33333  33333  333 


s  -  coordinates 

51  =  rj 

52  =  rj 

53  =  '*2 
H  “  ''6 

^5  "  '*4 
^6  '  '*4 

S7  - 


Weights 

Wj  *  0.12593  91805  448 
W2  = 

W3  =  Wj 

W4  «  0.13239  41527  885 

W5  =  W4 
Wg  =  W4 

=  0*225 


p1dC6in6nt  fields  synmetric  dbout  the  x  end  y  dxesa  Symnetry  dryuments 
allow  the  analysis  to  be  limited  to  one-quater  of  the  plate  and  hence 
reduce  the  deyrees  of  freedom  of  the  model.  Solviny  fracture  problems 
that  do  not  allow  the  use  of  symmetry  aryuments  requires  significantly 
laryer  amounts  of  computer  storage  or  more  sophisticated  storage  and 
equation  solving  routines  than  those  used  in  this  Investigation. 

3.3  Verification  of  the  Finite  Element  Model 

Distributions  of  the  stress  component'-;  near  the  crack  tip  for  an 
Infinite  center  cracked  plate  subjected  to  a  uniaxial  loading 
of  equal  to  1.0  ksi  and  a  crack  orientation  of  o  «'  0®  are  shown  in 
Flys.  3,5  to  3.13.  In  these  figures,  the  solid  lines  represent  results 
from  the  elasticity  solution  and  the  finite  element  results  for  yauss 
point  1  (Fig.  3.2)  are  represented  by  circles.  This  point  was  chosen 
since  it  is  the  nearest  integration  station  to  the  crack  tip  and  it  lies 
on  a  radial  line  through  the  centroid  of  the  element.  Similar  plots  for 
the  stress  distributions  at  other  gauss  points  are  presented  in  Appendix 
C. 

Figures  3.5  to  3.7  are  for  a  steel  plate.  For  this  isotropic 
problem,  the  finite  element  stress  results  compare  very  well  to  the 
elasticity  solution.  There  is  excellent  correlation  both  in  magnitude 
of  the  stresses  and  distribution  as  a  function  of  i). 

Figures  3,8  to  3.10  correspond  to  a  plate  manufactured  from 
Hercules  AS4/3501-6  graphite/epoxy  (material  properties  are  given  in 
Appendix  E).  For  this  problem,  the  fiber  orientation  angle,  e,  is  equal 
to  zero  degrees.  Overall,  there  is  correlation  in  the  stress 
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Steel  Plate 

a*o® 

Gauss  Point  I  0.004 
Oy  «I.O  KSI 


Fig.  3.5  OJstributlon  of  o  Near  the  Crack  Tip  In  a  Center  Cracked 
Steel  Plate* 


(KSI) 


0  (Degrees) 


Fig.  3.6  Distribution  of  o  Kear  the  Crack  Tip  In  a  Center  Cracked 
Steel  Plate.  ^ 


(KSI) 


Fig.  3 


7  Distribution  of  Near  the  Crack  Tip  In  a  Center 
Steel  Plate.  ^ 
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distribution  as  a  function  of  t  and  in  the  stress  magnitudes.  However, 
the  finite  element  prediction  of  the  magnitude  of  a^,  shown  in  Fig.  3.8, 
does  not  compare  as  well  as  the  other  stress  components.  The  lack  of 
better  agreement  for  this  stress  component  is  due  to  the  high  gradient 
that  exists  as  a  function  of  If  the  t/a  ratio  were  reduced  and  the 
concentration  of  elements  increased,  the  component  predicted  by  the 
finite  element  model  should  improve.  Since  the  gradients  in  th-f  stress 
components  for  graphite/epoxy  are  greater  than  that  of  steel,  there  is 
no  reason  to  expect  that  a  mesh  that  has  converged  for  an  isotropic 
fracture  problem  has  also  converged  for  a  highly  orthotropic  one. 

Similar  results  are  seen  in  Figs,  3.11  to  3.13  which  correspond  to 
a  plate  manufactured  from  AS4/3501-6  graphite/epoxy  with  a  fibers 
orientation  angle  of  90°,  The  distribution  in  stress  components  as  a 
function  of  (j.  compares  well  with  the  elasticity  solution  while  there  is 
some  discrepancy  in  the  magnitude  of  a^.  The  difference  in  magnitude  is 
less  than  that  of  the  o  »  0°  problem.  This  fact  is  not  surprising  since 
the  gradient  of  as  a  function  of  <|)  is  an  order  of  magnitude  smaller 
for  the  90°  problem. 

These  results  indicate  that  singular  finite  elements  can  be  applied 
to  anisotropic  fracture  problems.  However,  prior  to  their  application, 
a  detailed  study  on  the  convergence  of  singular  isoparametric  elements 
in  anisotropic  fracture  is  required.  Since  stress  gradients  in  aniso¬ 
tropic  materials  can  be  significantly  larger  than  those  in  isotropic 
materials,  the  mesh  building  guidelines  of  Freese  and  Tracey  [17]  are  no 
longer  applicable.  An  example  of  the  importance  of  the  1/a  ratio  on  the 
prediction  of  crack  tip  stresses  in  anisotropic  materials  is  shown  in 
Appendix  F. 
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AS4/350I-6 

a»0®  0*0* 

Gauss  Point  I  ’’^*0.004 
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Fig.  3.8  Distribution  of  Near  the  Crack  Tip  In  a  Center  Cracked 
Graphite/Epoxy  Plate  (e  -  0*). 
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Fig.  3.10  Distribution  of  Near  the  Crack  Tip  In  a  Center  Cracked 
Graphite/Epoxy  Plate  (e  -  0*). 
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AS4/350I-6  OTyf 


Fig.  3.11  Distribution  of  Near  the  Crack  Tip  in  a  Center  Cracked 
Graph1te/Epoxy  Plate  (e  =»  90®). 


(KSI) 


Fig.  3. 


Distribution  of  Near  the  Crack  Tip  in  a  Center  Cracked 
Graphite/Epoxy  Plate  (e  »  90*). 


Chapter  4 

OFF-AXIS  AND  UNIDIRECTIONAL  TENSILE  COUPONS 
In  this  section  analytical  results  are  presented  for  graphite/epoxy 
tensile  coupons  containing  center  cracks.  The  emphasis  of  the  analysis 
was  to  verify  the  numerical  results  of  Buczek  and  Herakovich  [5]  and  to 
develop  analytical  results  which  when  compared  to  existing  experimental 
results  would  test  the  credibility  of  the  Normal  Stress  Ratio  as  a 
viable  crack  extension  direction  criterion.  All  results  were  generated 
using  the  anisotropic  elasticity  solution  of  Section  (3.1).  Additional 
results  for  analyses  presented  in  this  chapter  can  be  found  in  Appendix 
D. 

4.1  Unidirectional  30®  Off-Axis  Tensile  Coupon 

Buczek  and  Herakovich  compared  several  crack  extension  direction 
theories  with  experimental  results  for  a  30®  off-axis  tensile 
coupon  (6  »  120®)  and  found  surprising  results.  Of  the  three  crack 
growth  direction  criteria  described  in  Section  (2.1),  only  the  Normal 
Stress  Ratio  agreed  with  the  experimental  value  of  300®.  Since  the 
finite  element  formulation  which  they  employed  was  not  developed  specif¬ 
ically  to  describe  the  stress  field  near  a  singularity,  it  was  desired 
to  test  the  theory  using  a  more  accurate  stress  analysis.  In  an  effort 
to  verify  their  results,  the  same  30®  off-axis  tensile  coupon  was  ana¬ 
lyzed  using  the  anisotropic  elasticity  solution. 

The  geometry  of  the  specimen  is  shown  in  Fig.  4.0.  The  crack  tip 
stresses  were  evaluated  for  an  applied  load  7^  »  1.0  ksi  at  an  r^  value 
of  0.005  inches.  An  r^  of  0.005  was  chosen  to  match  that  selected  by 
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Buczek  and  Herakovich.  The  distributions  of  the  crack  growth  direction 
criteria  as  a  function  of  t  are  shown  in  Figs.  4.1  through  4.3.  A 
comparison  of  the  predicted  crack  growth  directions  using  anisotropic 
elasticity  with  those  predicted  by  finite  elements  [5]  is  given  in  Table 

4.0. 

The  direction  of  crack  extension  predicted  by  the  Normal  Stress 
Ratio  is  very  consistent  in  both  the  analytical  and  finite  element 
solutions,  while  the  direction  of  crack  extension  predicted  by  the 
Tensor  Polynomial  and  Strain  Energy  Density  is  not  consistent.  The 
discrepancy  in  the  crack  growth  direction  predicted  by  the  Tensor 
Polynomial  and  Strain  Energy  Density  in  the  two  analyses  is  due  to  the 
finite  element  formulation  used  by  Buczek  and  Herakovich  [5].  Rather 
than  concentrating  on  approximating  the  stress  field  around  tne  crack 
tip,  their  formulation  emphasizes  evaluation  of  the  strain  energy 
released  during  crack  extension. 

It  is  important  to  note  that  the  conclusions  drawn  by  Buczek  and 
Herakovich  concerning  the  inconsistencies  of  the  Tensor  Polynomial  and 
Strain  Energy  Density  Criterion  in  predicting  the  direction  of  crack 
extension  are  supported  by  the  anisotropic  elasticity  solution.  The 
difference  in  the  value  of  predicted  by  the  two  methods  is  due  to  the 
improved  accuracy  of  the  anisotropic  elasticity  solution  relative  to  the 
conventional  finite  element  formulation. 

4.2  Influence  of  r^  on  the  Predicted  Direction  of  Crack  Growth 

Investigation  as  to  the  influence  of  r^  on  the  crack  growth 
direction  criteria  yields  interesting  results.  The  direction  of  crack 
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Fig.  4.1  Normal  Stress 
Coupon. 
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Fig.  4.3  Strain  Energy  Density  vs.  ♦  for  30*  Off-Axis  Graphite/Epoxy 
Coupon. 


Comparison  of  Anisotropic  Elasticity  to  Finite  Element  Results  for 
Direction  of  Crack  Extension  in  30*  Off-Axis  Tensile  Coupons 
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extension  predicted  by  the  Nornidl  Stress  Ratio  and  Strain  Energy  Density 
Criteria  Is  Independent  of  the  value  of  rg,  but  this  Is  not  the  case  for 
the  Tensor  Polynomial.  These  facts  are  shown  dramatically  In  Fig, 

4,4,  From  this  figure,  It  Is  apparent  that  a  small  change  In  the  value 
of  rg  can  affect  the  value  of  <p^  as  predicted  by  the  Tensor  Polynomial 
by  as  much  as  40®.  As  seen  In  Fig.  4.5,  the  change  In  the  distribution 
of  the  Tensor  Polynomial  as  a  function  of  ^  for  an  oroer  of  magnitude 
change  In  r^  Is  significant.  For  an  r^  of  5.0E-03",  the  Tensor 
Polynomial  predicts  the  direction  of  crack  growth  to  be  248®,  while  for 
an  rg  of  5.0E-04",  the  predicted  direction  Is  253®, 

4.3  Unidirectional  15®  Off-Axis  Coupons 

In  an  attempt  to  further  verify  the  consistency  of  the  Normal 
Stress  Ratio  Criterion,  the  experimental  configurations  of  Herakovich  et 
a1,  [26]  were  also  analyzed.  The  geometries  of  the  specimens  analyzed 
are  shown  In  Fig.  4.6.  The  experimental  observations  for  the  direction 
of  crack  extension  were  made  with  AS4/3501-6  graphite/epoxy.  The  ana¬ 
lysis  consisted  of  comparing  the  predicted  direction  of  crack  extension 
In  an  Infinite  plate  (e  »  105®)  with  an  applied  load  of  o’  ■  1,0  ksl  and 
an  rg  of  U.UU2  Inches  with  the  experimentally  observed  crack  extension 
direction.  An  r^  of  0,002  Inches  corresponds  with  the  location  of  gauss 
point  1  In  the  singular  Isoparametric  finite  element  solution. 

The  first  specimen  analyzed  has  a  crack  orientation  angle,  o,  of 
0®,  The  direction  of  crack  growth  predicted  by  the  Normal  Stress  Ratio 
for  this  sample  Is  286®.  This  corresponds  very  well  with  the  experimen¬ 
tal  value  of  285®  (parallel  with  the  fibers)  observed  by  Herakovich  et 
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F1g.  4.5  Influence  of  on  the  Distribution  of  the  Tensor  Polynomial 
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a1.  [26].  The  distribution  of  the  Normal  Stress  Ratio  as  a  function 
of  for  this  problem  is  shown  in  Fig.  4.7, 

For  the  second  specimen,  the  crack  was  aligned  perpendicular  to  the 
fibers  (o  =  15®).  The  distribution  of  the  Normal  Stress  Ratio  as  a 
function  of  for  this  problem  is  shown  in  Fig.  4.8.  The  Normal  Stress 
Ratio  predicts  the  direction  of  crack  extension  to  be  271®,  while  the 
experimental  results  of  Herakovich  et  al .  [26]  showed  crack  extension  to 
occur  at  270®.  Again  there  is  good  correlation  between  the  analytical 
model  and  the  experimental  results. 

Analysis  of  Figs,  4.7  and  4.8  yields  interesting  results.  For  both 
problems,  the  theoretical  prediction  of  differs  slightly  from  the 
experimentally  observed  value.  However,  for  both  cases  there  is  a 
strong  peak  in  the  distribution  of  the  Normal  Stress  Ratio  as  a  function 
of  ^  in  the  actual  direction  of  crack  growth.  This  fact  is  very  impor¬ 
tant,  The  Normal  Stress  Ratio  may  not  have  the  accuracy  to  predict  the 
direction  of  crack  excension  correctly  to  within  one  degree.  However, 
when  observed  graphically,  the  Normal  Stress  Ratio  represents  the  direc 
tion  of  crack  extension  exceptionally  well. 

4.4  Unidirectional  Lamina  Subjected  to  Mixed-Mode  Loadings 

In  order  to  study  crack  growth  in  fibrous  composites  under  more 
general  loading  conditions,  and  to  assist  the  experimentalist  in  verifi- 
cation  of  the  Normal  Stress  Ratio  Criterion,  unidirectional  lamina  were 
analyzed  for  various  mixed-mode  loadings  and  the  predicted  values 

of  determined.  Studies  of  thi s  nature  can  be  very  beneficial .  In 
addition  to  providing  theoretical  predictions  for  the  direction  of  crack 
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Fig.  4.7  Honnal  Stress  Ratio  vs.  ^  for  15*  Off-Axis  Graphite/tpoxy 
Coupon  (a  «  0*). 
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extension  under  mixed-mode  loading,  they  serve  as  a  guide  to  the 
experimentalist.  With  knowledge  of  the  predicted  response  In  hand,  one 
can  then  choose  the  most  appropriate  laboratory  experiments  necessary  to 
critically  test  the  hypothesis,  without  running  a  plethora  of  tests. 

The  geometry  and  loading  of  the  specimens  analyzed  Is  shown  In  Fig. 
4.9.  The  fiber  orientation  angle,  e.  Is  0®,  the  crack  orienta¬ 
tion,  o.  Is  0®  and  r^  was  chosen  to  be  0.002  Inches.  Mixed-mode  analy¬ 
sis  requires  the  definition  of  a  parameter,  X,  which,  Is  the  ratio 
of  F  to  T  .  Mode  I  loading  corresponds  to  X  »  »,  while  mode  II  Is 

y 

defined  as  X  =  0.  Figure  4.10  shows  the  predicted  direction  of  crack 
growth,  for  pure  mode  II  and  various  degrees  of  mixed-mode  loading. 

By  analyzing  Fig.  4.10  a  great  deal  of  Information  can  be  gained 
about  the  mixed  mode  response  of  0®  graphite/epoxy  material.  The  direc¬ 
tion  of  crack  extension  Is  strongly  dependent  on  x  for  small  value 
of  X  (l.e.,  <  2.0).  For  X  >  2.0  the  predicted  direction  of  crack  growth 
is  Independent  of  X.  For  pure  mode  II  loading,  the  predicted  crack 
extension  direction  Is  348®.  In  contrast,  a  mixed-mode  loading  of  X  »  1 
generates  crack  extension  at  358®,  while  pure  mode  I  generates  crack 
growth  at  0®.  From  these  results,  an  Ideal  laboratory  test  to  verify 
the  Normal  Stress  Ratio  would  be  pure  mode  11  loading.  The  predicted 
direction  of  crack  extension  for  this  loading,  348®,  clearly  Implies 
fiber  breakage.  A  very  poor  test  would  be  any  mixed-mode  test  that  was 
not  predominantly  mode  II  because  the  direction  of  crack  extension 
appears  to  be  dominated  by  the  mode  I  response  when  X  >  2.0. 


Chapter  5 

CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  has  been  concerned  with  the  development  of  models  to 
predict  the  crack  tip  stress  field  In  unidirectional  composite  materi¬ 
als,  Investigation  of  various  crack  yrowth  direction  criteria,  and  the 
effect  of  biaxial  loading  on  the  direction  cf  crack  extension. 

For  the  cases  analyzed,  the  Normal  Stress  Ratio  has  proven  to  be  a 
consistent  crack  extension  direction  criteria.  Two  advantageous  charac¬ 
teristics  of  the  Normal  Stress  Ratio  are,  the  presence  of  a  strong  peak 
In  the  direction  of  crack  extension  when  plotted  as  a  function  of  t  and 
the  Independence  of  on  the  choice  of  r^.  Though  the  Normal  Stress 
Ratio  may  not  have  the  accuracy  to  predict  the  direction  of  crack  exten¬ 
sion  to  within  one  degree,  when  observed  graphically,  this  criterion 
predicts  the  direction  of  crack  extension  exceptionally  well.  The 
Independence  of  on  the  choice  of  r^  allojvs  the  prediction  of  crack 
extension  direction  in  unidirectional  composites  without  prior  knowledge 
of  a  proper  r^  value. 

This  study  has  confirmed  the  discrepancies  found  by  Buczek  and 
Herakovich  [5]  regarding  the  Tensor  Polynomial  and  Strain  Energy  Density 
crack  growth  criteria.  The  inability  of  these  theories  to  predict  the 
correct  direction  of  crack  growth  was  not  due  to  an  inaccurate  finite 
element  stress  analysis,  but  due  to  limitations  In  the  theories. 

Application  of  the  anisotropic  elasticity  solution  to  the  analysis 
of  unidirectional  composites  containing  cracks  has  been  very  advantage¬ 
ous.  The  speed  and  reduction  in  cost  of  analyses  based  on  the  aniso¬ 
tropic  elasticity  solution  relative  to  finite  element  solutions  is 
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phenomenal.  This  access  to  an  Inexpensive  analytical  solution  allows 
the  performance  of  detailed  parametric  studies.  Through  these  paramet¬ 
ric  studies,  critical  tests  can  be  determined  to  better  verify  and 
evaluate  crack  direction  criteria. 

It  has  been  shown  that  singular  Isoparametric  finite  elements  can 
also  be  used  to  approximate  crack  tip  stresses  In  anisotropic  materi¬ 
als,  Though  the  problems  treated  in  this  study  where  heuristic  In 
nature,  they  have  proven  that  singular  Isoparametric  finite  elements  can 
be  used  to  solve  anisotropic  fracture  problems.  Finite  element  solu¬ 
tions  are  particularly  valuable  for  problems  that  cannot  be  addressed  by 
the  anisotropic  elasticity  solution,  1,e,,  problems  with  finite  size  or 
complex  geometries. 

Though  this  study  has  proven  the  potential  usefulness  of  singular 
Isoparametric  finite  elements  In  anisotropic  fracture,  additional 
research  Is  required.  Particular  emphasis  should  be  given  to  the  devel¬ 
opment  of  guidelines  which  will  assist  the  analyst  In  mesh  building, 
i.e,,  proper  t/a  ratios  and  element  concentrations. 

Additional  experimental  research  is  required  to  further  verify  the 
Normal  Stress  Ratio  Criterion.  Particular  consideration  should  be  given 
to  performing  tests  under  mode  II  loading.  The  Normal  Stress  Ratio 
predicts  fiber  breakage  for  pure  mode  II  loading  of  a  unidirectional 
graphite/epoxy  composite  with  fiber  orientation  and  crack  orientation  of 
0*.  Correspondence  between  experimentally  observed  and  analytically 
predicted  values  of  for  tests  Involving  fiber  breakage  would  add 
significant  credibility  to  the  Nonnal  Stress  Ratio  Criterion. 


Lastly,  to  fully  verify  the  Normal  Stress  Ratio  Criterion  other 
material  systems  must  be  studied.  So  far,  all  applications  of  the 
Normal  Stress  Ratio  Criterion  have  been  with  graphite/epoxy.  Material 
systems  representing  various  degrees  of  anisotropy  from  Isotropic  to 
graphite/epoxy  need  to  be  Investigated. 
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APPENDIX  A 


1-2  Coordinate  System 

f{o^)  »  FjOj  +  F2a2  +  F^^Oj^  +  F2202^  +  2FJ2O1O2  +  Fg^a^^ 

Mathematical  Expressions  for  the  Strenyth  Tensor  Terms  as  a 
Function  of  Measured  Strengths 
Fi  •  {\/Sj  *  l/Xj.) 

Fg  -  {1/Yx  + 

hi  •  -i/(XtXc) 

F22  -  -i/{Yt^c) 

F66  • 


X-Y  Coordinate  System 

f(ai)  •  F^o^  +  FyO^  ♦  FjT^^  +  *  Vy^  *  '' 

■*■  ^^xy®x°y  '*■  ^'^xS°x^xy  ^^yS®y^xy 
Transformation  of  Strenyth  Tensor  Terms 


SS^xy 


2 


Fj^  »  m^Fj  +  n^F2 
Fy  ■  n^Fj  »  m2F2 
Fj  *  2mnFj  -  2fnnF2 

Fxx  -  +  "^^22  2m2n2Fi2  +  m^n^Fgg) 

Fyy  -  (n^Fjj  +  m^F22  +  *  ">^"^^66  ) 

Fgs  ■  (Am^n^Fj^  +  4ffl2n^F22  - 

Fxy  •  (m^n^Fjj  +  m^ii^F22  +  (m^  +  n^)  Fj2  - 

Fxs  ■  (ZM^nFjx  -  2mn3F22  +  2(irn3  ,  n,3n)Fj2  +  (mn^  -  m^n}Fgg) 

Fys  “  (2mn^F^j  -  2m^nF22  +  2(m^n  -  mn^)Fj2  +  (m^n  -  mn^)Fgg) 


where  m  *  cos  6  and  n  »  sin  6. 
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appendix  U 


To  provide  a  better  understanding  of  singular  Isoparametric  finite 
elements,  the  form  of  the  singularity  along  the  line  1-2  (s  -  0)  In  Fig 
8.1  will  be  calculated. 

For  a  six  noded  Isoparametric  triangular  element,  the  following 
relations  hold: 

"■  j,  (1) 

“  ■  (2) 

where 

h^  Is  the  shape  function  for  the  1^^  node 
Is  the  x-coordinate  of  the  1^^  node 

u^  Is  the  x-directlon  displacement  of  the  1^^  node 
X  Is  the  x-coordInate  mapping  function 
u  Is  the  x-directlon  displacement  Interpolation  function 

Evaluation  of  the  shape  functions,  along  the  line  1-2  (s  •  0), 
considering  nodes  1,  2  and  4  only  yields: 

hj  »  1  -  3r  +  2r^ 
hg  “  -  r  +  2r 
^4  “  4r  -  4r^ 

Substitution  of  these  shape  functions  Into  Eq.  1  yields: 
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Conventional  Triangular  Element 


4  2 


Fig.  B.l  Node  Locations  <n  Conventional  and  Singular  Isoparametric 
Triangular  Elements. 
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(9) 
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Analysis  of  Eq,  9  reveals  that  the  strain  singularity  along  the  line  1-2 

Is  of  the  order  —  ,  which  Is  the  required  singularity  for  fracture 
/r 

mechanics  applications. 


APPENDIX  C 

The  followiny  figures  show  the  distribution  In  the  stress 
components  as  a  function  of  ^  for  the  seven  gauss  points*  This 
Information  Is  Included  to  give  a  better  understanding  of  the  stress 
field  distribution  through  the  element.  The  problem  under  consideration 
Is  a  center  cracked  AS4/3501-6  plate  (6  >  90”,  2a  »  0.5")  subjected  to  a 
far  field  stress  of  1.0  KSl. 


(KSI) 
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AS4/350I-6 
a«0®  0*90® 

6au*e  Point  I  0.004 
Oy  » 1.0  KSI 


Fig.  c.l 


♦  at  Gauss  Point  I  for  a  Center  Cracked  Graphite/Epoxy 


(KSI) 


AS4/350I-6 
a«0®  0*90® 

Cause  Point  I  0.004 
Oj  •  1.0  KSI 


Fig.  C.2  Oy  vs.  «  at  Gauss  Point  1  for  a  Center  Cracked  Graphite/Epoxy 
Plate. 


(KSI) 


(KSI) 
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Fig.  C.5  o„  vs.  4  at  Gauss  Point  2  for  a  Center  Cracked  Graphite/Epoxy 


(KSI) 


AS4/350I-6 
a«0*  0«9O* 
Gausa  Point  2 
aj  •  1.0  KSI 


0  (Dagreas) 


Txy  vs.  4  at  Gauss  Point  2  for  a  Center  Cracked  Graphite/Epoxy 


0  (Degrees) 


Fig.  C.7  vs.  ^  at  Gauss  Point  3  for  a  Center  Cracked  Graphite/Epoxy 
Plate. 


Fig.  C.8  Oy  vs.  ♦  at  Gauss  Point  3  for  a  Center  Cracked  Graphite/Epoxy 
Plate. 


(KSI) 


0  (Degrees) 


Fig.  C.ll 


Oy  v$.  ♦  at  Gauss  Point  4  for  a  Center  Cracked  Graph ite/Epoxy 


-  —  Anisotropic  Elosticity  Solution 


0  30  60  90  120  160  180 

0  (Degrees) 

Fig.  C.12  vs.  ^  at  Gauss  Point  4  for  a  Center  Cracked 
Graphite/Epoxy  Plate. 
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Fig.  C. 
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AS4/350I-6 


13  Ojj  vs.  ^  at  Gauss  Point  5  for  a  Center  Cracked  Graphite/Epoxy 
Plate. 


AS4/350I-6 
a«0»  d»90* 

6aus«  Point  5  ^^>0.09 
^  -I.OKSI 


— Anisotropic  Elasticity  Solution 
"  o  Finite  Element 


0  30  60  90  120  160 

0  (Degrees) 


Fig.  C.14  Oy  vs.  «  at  Gauss  Point  5  for  a  Center  Cracked  Graphite/Epoxy 
Plate. 
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AS4/350I-6 

a»o*  0*90* 

Goum  Point  5  %b0.09 
^  >1.0  KSI 


o  o 


_ — Anisotropic  Elasticity  Solution 
o  Finite  Element 


-1.0 


0  30  60  90  120  160  ISO 

0  (Degrees) 


Fig.  C.15  Tjjy  vs.  ♦  at  Gauss  Poirtt  5  for  a  Center  Cracked 
Graphite/Epoxy  Plate. 
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Anisotropic  Elasticity  Solution 
o  FInit®  Element 
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AS4/350I-6 
a*0®  0-9O* 

Gcusa  Point  6  '^■0.03 
Ol  ■1.0  KSI 


—  Anisotropic  Elasticity  Solution 
o  Finite  Element 


O  30  60  90  120  160  180 

0  (Degrees) 


Fig.  C.17  o  vs.  «  at  Gauss  Point  6  for  a  Center  Cracked  Graphite/Epoxy 
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Flg^  C.19 


vs.  *  at  Gauss  Point  7  for  a  Center  Cracked  Graphite/Epoxy 


/xod3/9)(i|d9ag  p9:|3SJ3  J9)ua3  «  uo^  i  )u^Od  ssnsQ  ♦ 


(S90j5eQ)  0 


Axtitotropic  Elasticity  Solution 
o  Rnite  Elomsnt 


APPENDIX  0 


This  section  is  included  to  pro-ide  additional  information 
concerning  the  prediction  of  crack  extension  In  unidirectional 
composites.  The  figures  presented,  provide  Insight  into  the  crack  tip 
stress  states  and  characteristics  of  the  crack  extension  direction 
criteria. 
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0  ( 0«9r9tt) 

Fig.  0.1  Influence  of  on  the  Distribution  of  In  a  Center  Cracked 
6raph1te/Epoxy  Plate  (a  »  120*), 
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(KSI) 
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T300/5208 
a«0*  8«I20® 

ro  ■  0.002  in 
5*y-l.0KSI 


0  (0«gre9t) 


Fig.  0.2  Influence  of  on  the  Distribution  of  o  In  a  Center  Cracked 
Graphite/Epoxy  Plate  (e  ■  120*), 


(Odgress) 


3  Influence  of  on  the  Distribution  of  t 
Graphite/Epoxy  Plate  (0  =  120*) . 


0  ( Dogrses) 


D.7 


Tensor  Polynomial  vs.  «  for  a  Center  Cracked  30 
Graphite/Epoxy  Plate  Under  Pure  Shear. 


0  ( DsQraes) 


Fig.  0.12  Tensor  Polynomial  vs.  a  for  a  Center  Cracked  15 
Graphite/Epoxy  Plate. 


0  ( Degrees) 


13  Strain  Energy  Density  vs.  4  for  a  Center  Cracked  15 
Graphite/Epoxy  Plate. 


0  (Degrees) 


vs.  «  for  a  Center  Cracked  Graphite/Epoxy  Plate  («  ■  0 


p  (Oogrees) 


16  Nonnal 
Plate 


Stress  Ratio  vs.  a  for  a  Center  Cracked  Graphite/Epoxy 
o  »  0*  and  e  «  0*). 


<P  { Degrees) 


17  Tensor  Polynomial  vs.  a  for  a  Center  Cracked  Graphite/Epoxy 
Plate  (o  ■  0*  and  e  ■  0*). 


Strain  Er 


g.  0.18  Strain  Energy  Density  vs.  e  for  a  Center  Cracked 
Graphite/Epoxy  Plate  (a  -  0‘  and  e  «  0*). 
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APPENDIX  E 


AS4/3501-6  Graphite/Epoxy 

Ej  ■  21.6  HSI 
Eg  •  1.96  M5,I 
Gjg  •  0.83  HSI 
Vj2  *  0.28 
Xy  -282  KSI 
Yj  -  10  KSI 
S  -  14.2  KSI 


T300/5208  Graphite/Epoxy 

El  •  10.2  MSI 
Eg  ■  1.56  MSI 
Gig  >  0.82  HSI 
^12  ■  0,24 
Xy  •  219.5  KSI 
Xq  •  -246  KSI 
Yy  »  6.35  KSI 
Yq  »  -23.8  KSI 
S  «  12.6  KSI 
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APPENDIX  F 


This  section  Is  Included  to  provide  additional  Information 
concerning  the  convergence  of  the  singular  finite  element  solution. 
Results  presented  are  for  the  finite  element  mesh  shown  In  Fig.  F.l. 
For  this  mesh  the  t/a  ratio  Is  l.O  and  the  2a/w  ratio  Is  0.1.  The 
problem  under  consideration  Is  a  center  cracked  AS4/3501-6  plate  (e  ■ 
90“  and  a  -  0“)  subjected  to  a  far  field  stress  7^  •  l.o  tCSI. 

Distributions  of  the  stress  components  near  the  crack  tip  as  a 
function  of  «  are  given  In  Figs.  F.2  to  F.4.  In  these  figures,  the 
solid  Vine  represents  results  from  the  anisotropic  elasticity  solution 
and  the  finite  element  results  for  gauss  point  1  are  represented  by 
circles.  Analysis  of  the  figures  Indicates  the  finite  element  solution 
adequately  approximates  the  distribution  of  the  stress  components  as  a 
function  of  However,  there  Is  a  large  discrepancy  In  the  magnitude 
of  the  stresses  predicted  by  the  finite  element  solution  and  those 
predicted  by  the  anisotropic  elasticity  solution. 

Comparison  of  these  results  with  the  results  of  Chapter  3  shows  the 
Importance  of  the  t/a  ratio  on  the  prediction  of  crack  tip  stresses. 

The  finite  element  mesh  used  In  Chapter  3  has  a  smaller  t/a  ratio  (1/a  * 
0.1)  and  hence  better  approximates  the  crack  tip  stress  field. 
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Plate. 


♦  at  Gauss  Point  1  for  a  Center  Cracked  Graphite/Epoxy 


Fig.  F.2 
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Gauss  Point  1  for  a  Center  Cracked 
ly  Plate. 


VIRGINIA  TECH  CENTER  FOR 
COMPOSITE  MATERIALS  AND  STRUCTURES 


The  Center  for  Composite  Materials  and 
Structures  is  a  coordinating  organization  for 
research  and  educational  activity  at  Virginia 
Tech.  The  Center  was  formed  in  1982  to 
encourage  and  promote  continued  advances  in 
composite  materials  and  composite  structures. 
Those  advances  will  be  made  from  the  base  of 
individual  accomplishments  of  the  thirty-four 
founding  members  who  represent  ten  different 
departments  in  two  colleges. 

The  Center  functions  by  means  of  an 
Administrative  Board  which  is  elected  yearly. 
The  general  purposes  of  the  Center  include: 

•  collection  and  dissemination  of  informa¬ 
tion  about  composites  activities  at  Virginia 
Tech, 

•  contact  point  for  other  organizations  and 
individuals, 

•  mechanism  for  collective  educational  and 
research  pursuits, 

•  forum  and  mechartism  for  internal  inter¬ 
actions  at  Virginia  tech. 

The  Center  for  Composite  Materials  and 
Structures  is  supported  by  a  vigorous  program 
of  activity  at  Virginia  Tech  that  has  developed 
since  1963.  Research  expenditures  for  investiga¬ 
tions  of  composite  materials  and  structures  total 
well  over  five  million  dollars  with  yearly 
expenditures  presently  approaching  two  million 
dollars. 


Research  is  conducted  in  a  wide  variety  of 
areas  including  design  and  analysis  of  compo¬ 
site  materials  and  composite  structures, 
chemistry  of  materials  and  surfaces,  characteri¬ 
zation  of  material  properties,  development  of 
new  material  systems,  and  relations  between 
damage  and  response  of  composites.  Extensive 
laboratories  are  available  for  mechanical 
testing,  nondestructive  testing  and  evaluation, 
stress  analysis,  polymer  synthesis  and  character¬ 
ization,  material  surface  characterization, 
component  fabrication  and  other  specialties. 

Educational  activities  include  eight  formal 
courses  offered  at  the  undergraduate  and 
graduate  levels  dealing  with  the  physics, 
chemistry,  mechanics,  and  design  of  composite 
materials  and  structures.  As  of  1982,  some  33 
Doctoral  and  37  Master's  students  have 
completed  graduate  programs  and  several 
hundred  Bachelor-level  students  have  been 
trained  in  various  aspects  of  composite 
materials  and  structures.  A  significant  number 
of  graduates  are  now  active  in  industry  and 
government. 

Various  Center  faculty  are  internationally 
recognized  for  their  leadership  in  composite 
materials  and  composite  structures  through 
books,  lectures,  workshops,  professional  society 
activities,  and  research  papers. 
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and  Operations  Research 

Raphael  T.  Haltka 

Hal  F.  Brinson 

|oel  K  Nachlas 
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Materiris  Engineering 

Eric  R.  Johnton 

Daniel  Frederick 

David  W.  Dwight 
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Robert  A.  Heller 

D.  P.  H.  Hasselman 

Donald  G.  Baird 

Edmund  G.  Henneke,  11 
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Carl  T.  Herakovkh 
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Qvii  Engineering 

Don  H.  Morrb 
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S.  W,  Zewari 
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Inquiries  should  be  directed  to: 

Center  for  Composite  Materials  &  Structures 
College  of  Engineering 
Virginia  Tech 
Blacksburg,  VA  24061 
Phone:  (703)  961-4969 
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